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Abstract. The work in this article is inspired by a classical problem: the statistical 
physical properties of a closed polymer loop that is wound around a rod. Historically 
the preserved topology of this system has been addressed through identification of 
similarities with magnetic systems. We treat the topological invariance in terms of a 
set of rules that describe all augmentations by additional arcs of some fundamental 
basic loop of a given winding number. These augmentations satisfy the Reidemeister 
move relevant for the polymer with respect to the rod. The topologically constrained 
polymer partition function is now constructed using the combinatorics of allowed arc 
additions and their appropriate statistical weights. We illustrate how, for winding 
number 1, we can formally derive expressions for lower and upper bounds on the 
partition function. Using the lower bound approximation we investigate a flexible 
polymer loop wound between two slits, calculating the force on the slit as well as the 
average numbers of arc types in dependence of slit width and separation. Results may 
be extended to higher winding numbers. The intuitive nature of this combinatoric 
scheme allows the development of a variety of approximations and generalisations. 
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Entanglements occur naturally in polymeric systems, and are related to topological 
constraints. Indeed, the fact that different polymer strands cannot pass through each 
other is manifest in two observations: a) there exists an excluded volume interaction in 
real polymers that leads to self-avoidance, and b) for closed loops topological states must 
be conserved. We will focus on the latter statement, and concern ourselves with the 
associated conhgurational constraints. Such constraints determine which conhgurations 
of the polymeric system are topologically equivalent to each other, and thus restrict the 
polymer conhgurations over which we must sum to calculate the partition function for 
a given topological state. Mathematically the notion of topological equivalence may be 
captured (at least partially) by topological invariants. The role of topological constraints 
in polymers remains an important issue in various systems (see reviews by Kholodenko 
and Vilgis (^, and, more recently, mil)- Entanglement of synthetic or biological 
macromolecule loops continues to be treated in analytical and computational modelling 
(for example, in [5||^) and is deemed to be particularly relevant for localisation of DNA 
in cells [^. Although computer simulations have driven results strongly, there is still a 
need for expanding the range of analytical tools to deal with entangled chains. 

As early as 1961, Frisch and Wasserman considered topological isomerism in 
chemical systems, investigating knotted and unknotted links, loops and rings that occur 
in certain chemical molecules. Soon thereafter (1967), Edwards 10 explicitly pointed 
out the importance of topological constraints to polymeric systems, and that such 
constraints need to be included into the statistical mechanics of polymers. In that article, 
the specihc example of a polymer wound around a rod is considered, and the winding 
number is identihed as a suitable topological invariant that categorises topologically 
distinct conhgurations. The winding number simply represents the number of times 
the polymer winds around the rod before closing on itself. Viewing this problem in a 
projection along the rod, one may capture the planar winding number (angle) of the 
polymer around a hxed point (i.e., the projection of the rod) through the integral 
xy — yx 


+ y^ 


ds, 


( 1 ) 


x^ 

where s is the arc-length parameter of the strand, and x = ^ etc. Naturally this 


winding number cannot be altered once the polymer loop has been wound around the 
rod and closed. In |[^ this constraint is included into the path integral partition function 
(probability distribution) of the polymer inside a delta function that is exponentiated 
through the introduction of auxiliary helds. The resulting action is related to that of 
a magnetic system, and is treated under certain approximations. In the same year as 
Edwards, the same physical system was investigated by Prager and Frisch |[^. This 
problem can also be addressed through the introduction of a tailored potential that 
approximates an interaction at short distances with the rod, and is included in the path 
integral through an additional Boltzmann factor 


12,13 . We shall, however, focus on 


topological invariance as the basis of our approach. 
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A further simple invariant is the linking number, which tells us how many times 
two distinct loops a and b are wound around each other. The linking number is given 
by the Gauss integral 

ra - n 


ra{Sa) X fb{Sb) 


nl 


dSfi dsb^ 


( 2 ) 


where the position vectors r and their derivatives r are parametrised by the arc-lengths 
Sa and Sb- This invariant may also be included as a constraint in the polymer path 
integral (see, for instance, [1, 10, 14-l6|). 

Winding and linking numbers (and other basic invariants) do not capture all 
topological constraints in polymer systems. More complex higher order invariants (which 
may represent more detailed topological information) exist and may, in principle, also 
be included as constraints in path integrals. Indeed, the field theories that thus arise 
from topologically constrained polymer systems have been studied extensively [1,14 16 


There is also a deep connection between polynomial knot invariants and quantum 
field theories. The reader is referred to the seminal work by Witten 17 where Jones 
polynomials are investigated in the setting of Yang-Mills theory. This work opened the 
door to subsequent extensions such as perturbative approaches - see, for instance, 18 


Typically, however, simple invariants such as winding numbers and linking numbers 
have been considered in the context of polymer path integrals since more complex 
invariants become mathematically tedious to handle. The conservation of these 
topological invariants has been shown to relate to symmetry transformations that 


ultimately manifest in local gauge invariance in such field theories 14,15 


In this article we shall address the often studied problem of a polymer wound 
around an inhnitely long obstacle in a plane. In doing so, we shall not consider any 
self-entanglements of the loop, but simply concern ourselves with the topology of the 
loop relative to the rod. Consequently the mathematical intricacies of higher order 
invariants will not be of bearing here: we need look no further than winding numbers to 
address this physical system. Perhaps it is (in part) for this reason that the “polymer 
wound around a rod” has been studied so extensively. More recently (in 2003) 
Grosberg and Frisch |19l| presented various modifications and extensions 


of Edwards’ original results, both in the quenched (constrained partition 
sum) and the annealed (probability distributions of winding angles) settings. 
These include confining the polymer-and-rod system to a cavity, and winding 
the polymer around a disc. Our aim in this article is similar that of and 11 : 
we wish to find the partition function of a polymer wound around an obstacle. Our 
partition function, obtained through a different calculational approach, is 
then used to study physical quantities for various geometries. 

We employ a strategy that differs significantly from the path integral 
schemes cited above, namely to evaluate invariant knots by a combinatoric 
scheme. This may be done in terms of enumerations of braids, as shown by 
Nechaev and co-workers |20| , |21| . We also develop a combinatoric scheme, 
and use this for enumerating configurations subject to a winding number 
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constraint in particular. In principle the sequence of allowed topological 
configurations may be generated braid-theoretically. Since we, however, wish 
to calculate the partition function of polymer degrees of freedom and related 
quantities, it further is necessary to couple these degrees of freedom in a more 
careful manner to the combinatorics. This is presented in two parts. Firstly 
we outline how configurations may be labelled according to piercings that the 
strand makes through a plane. Then we consistently approximate the sum of 
all unique but topologically equivalent configurations comprised of connected 
arcs restricted to half-space. This method avoids the complications brought 
about by delta-function constraints which would be necessary to enforce the 
winding number in various alternate scenarios. This enables us to find upper 
and lower bounds for the free energy. In principle, this formalism allows for 
the description of various types of polymer chains, e.g., Gaussian or semi- 
fiexible. Secondly, we illustrate the calculation of the partition function and related 
averages for various winding scenarios, e.g., winding between two slits. We are able to 
calculate forces and ratios of arc types for confining geometries. To this end we consider 
Gaussian chains and their associated probability distributions, and calculate statistical 
quantities of interest. 


2. Winding a polymer around a rod 


The basis of our problem is a topological obstacle around which a polymer strand is 
wound. We start by presenting some simple examples of loops wound around a rod, 

and then illustrate different configurational modifications / augmentations 
that do not alter the topology of these basic loops. Such procedures are 
essentially braid manipulations on two strands which may be represented 
in terms of braid groups — see, for instance [22[ , |23| . Since we immediately 
couple the configurations of this quenched scenario to polymer degrees of 
freedom, however, the braid group relations alone are not sufficient for 
our enumeration procedure. This is discussed in detail in [Appendix A 
reference to 


m 


20,21 


2.1. Example of the basic loop, winding number w = 1 

Consider an infinitely long rod that is placed along the y axis of a system of axes in 
Suppose now that an open polymer strand is wound around this rod and then closed on 
itself to form a closed loop. The first natural question to ask is how often the polymer 
is wound around the rod. Indeed, this number distinguishes topologically distinct 
configurations of the polymer, and is appropriately known as the winding number, w. 
For the remainder of the article, the scenario in Figure with the minimal number of 
arcs will be referred to as the basic loop. If we take this configuration with w = 1 (the 
simplest case where strand is wrapped around the rod only once) we cannot deform 
or alter this configuration to obtain one where w ^ 1 without physically breaking the 
polymer strand. In this sense the winding number is a topological invariant of the 
particular configuration created when closing the open strand after w windings. 
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Figure 1: Closed polymer loop, w = 1. 


Let us divide the complete polymer loop into sub-arcs, the division occurring 
whenever the a;j/-plane is pierced by the polymer. For the example in Figure the 
entire polymer loop may be viewed as consisting of a “sequence” of two polymer arcs, 
each constrained by the x|/-plane to a half-space w.r.t. the positive / negative z axis. 
In this hgure we have labelled each of these arc-segments with a T whose subscript is 
H— if the arc begins in the a; > 0 half-plane {{x,y) : x > 0, y E M} and ends in the 
a; < 0 half-plane {{x,y) : x < 0, y E M}. The subscript —h applies to the opposite case. 
The Ts themselves will later represent the probability distributions for the half-space 
restricted sub-arcs. They are functions of the (planar) beginning and end co-ordinates 
of the respective polymer segments and of the arc-lengths of the segments. For now we 
will simply use this notation to represent sequences of such sub-arcs, and demonstrate 
how we may capture the topology of any polymer loop as a composition / sequence of 
such Ts. In this spirit, we represent the simple closed loop of Figure [T] symbolically by 
the sequence T+_T_+ (or alternatively by the cyclic permutation T__|_T+_). Essentially 
the sequence of subscripts indicates how one would follow the polymer strand around 
the rod from one piercing of the plane to another. An orientation convention (indicated 
by the arrows in Figure is chosen without loss of generality. The partition function 
corresponding to Figure in less compact notation would be 

Z = [ T+_(ri,f 2 )T_+(r 2 ,ri) = [ T_+(r 2 , n) T+_(ri, r 2 ), (3) 

Jv Jv 

where the position vectors label the piercings of the xy-plane and integration is over the 
relevant domain 

V = {xiE (0, oo); X2E (0, -cx)); yi, y2 E {-00, cx))}. (4) 

Each T{fi,fj) represents the statistical weight of a polymer arc restricted to half-space 
with appropriate initial and hnal positions in the plane. T depends on the nature of 
the specihc polymer. In principle the methods shown here can be applied to Gaussian 

as well as semiflexible polymers, etc. The symbolic sequence _T__(_ thus represents 

the integrand of the partition function for the basic loop. This notion will be clarihed 
in section We shall now extend this picture (and the symbolic notation) to higher 
winding numbers. 
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Figure 2 : Closed polymer loop, w = 2. 


Again the complete polymer strand may be decomposed into a sequence of conhned 
sub-arcs, each living in either the ^ > 0 or the ^ < 0 half-space. Following the 
notation above, the case w = 2 can thus be described symbolically by the sequence 
r+_T_+T+ _r__|_ (or any cyclic permutation thereof, depending on the choice of reference 
point for labeling). The partition function for Figure]^ would be 

Z = [ F+_(ri, Fa) T_+(r 2 , Fg) r+_(F 3 , F 4 ) T_+(F 4 , Fi) 

Jv 

= f T_+(F 2 ,F 3 )T+_(F 3 ,F 4 )T_+(F 4 ,Fi)T+_(Fi,F 2 ) etc. (5) 

Jv 

Clearly Z is invariant under cyclic permutation of the factors in the integrand in equation 
(|^. Consequently one could just as well label Figure with any cyclic permutation of 
T+_T_+T+_T_+. The order of the Ts does, however, matter, since the arguments of 
consecutive Ts (i.e., Fj, F+i etc.) must match up. This may be viewed in analogy to 
operator multiplication. 

The strategy in section will be to “diagonalise” the integral above, so that we 
may write symbolically T+_T_+T+_T_+ = (T+_T__|_)^. Analogously loops wound w 
times are expressed as (T+_T_+)"' in this compact notation. 

The examples considered thus far only show limited conhgurations associated with 
specific winding numbers, since they are composed of sub-arcs that cross from one side 
of the rod to the other. Other permissible configurations (that maintain the winding 
number) can include sub-arcs that remain on one side of the rod. In the next section 
we illustrate how simple loops such as those in Figures [T] and may be augmented in 
this way. Some combinatoric rules will be established on the symbolic level of sequences 
of Ts. The connection of these combinatoric sequences to a complete partition function 
will be made in section We shall now focus on the case of w = 1, since the partition 
functions for higher winding numbers are generated from powers of the basic loop. 
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2.3. Augmenting the basic loop: insertion of sub-arcs 

Let us return to the basic loop from Figure [T] with w = 1. We note that it is possible 
to augment or “decorate” this simple loop with more half-space constrained sub-arc 
segments. This process is subject to a Reidemeister move of the second type (see Figure 
of the polymer strand relative to the rod, viewed in a side-on projection along the x 
axis. 



Figure 3: Type two Reidemeister move of the polymer (thin) relative to the rod (thick). 


Since the Reidemeister moves do not alter a particular topological state for closed 


strands 24 , this augmentation does not alter the winding number but simply introduces 


more piercings of the xy-plane. Clearly the number of Ts equals the number of piercings 
for a particular conhguration. Introducing additional piercings / arc-segments may be 
done in two ways: 

(i) We may insert two more sub-arcs, each beginning and ending in the half-plane 
{{x,y) : X > 0, y E M}, one living in half-space z < 0 and the other in 2 ; > 0: 



Figure 4: Closed polymer loop, w = 1, additional constrained arc-segments. 


We describe the sequence of arcs in Figure with the symbolic sequence 
T+_T_+r++r+_|_ (or any cyclic permutation). Here the subscript indicates a 
sub-arc beginning and ending in the half-plane {{x,y) : x > 0, y E M}. 

(ii) Given Figure 0 we observe that we could also take one of the T++ sub-strands and 
“pull it acros^ the rod, as follows: 
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Figure 5: Closed polymer loop, w = 1, a further augmentation. 

The diagram shown here could be described by the sequence T+_T_+T+_|_T+_T_T_+ 

(or any cyclic permutation thereof). 

These two augmentation procedures form the basis for a set of combinatoric rules that 
govern what sequences are derivable from the basic loop with w = 1. Let us return 
once more to said basic loop in Figure From the examples above it is clear that this 
is the simplest loop for two reasons: (i) it has the smallest possible winding number 
w = 1, and (ii) it has the smallest number of piercings of the xy-plane. Clearly the 
two examples in Figures and are topologically equivalent to that in Figure they 
have the same winding number. However, they have more piercings of the a:|/-plane 
because more sub-arcs were inserted. The inclusion of these augmenting sub-arcs was 
topologically consistent: the winding number was not altered and the strand was not 
broken. We continue now by stating concretely what rules govern the augmentation of 
basic loops through insertion of sub-arcs in such a way that the basic topology (i.e., 
their winding numbers) is conserved. 

2.4- Condensed notation 

The statistical weights of the sub-arcs considered here are symmetric around z = 0. It 
is, however, important to distinguish between Ts that cross the rod and those whose 
two piercings of the plane are on the same side of the rod. We introduce the following 


shorthand. 




1 

+ 

(6) 

and 




T++, T — —)■ Ts, 

(7) 


where the subscripts c and s refer to “crossing” and “same side”, respectively. For 
Figure]^ for example, we may write 
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Of course it is implicit that a string be uninterrupted in its subscripts: two consecutive 
Ts of the form 

with a,/?, 7,5 e {+,-} (9) 

must be such that /? = 7 . If this were not the case one would have a broken strand since 
consecutive sub-arcs in different half-planes cannot be connected due to the rod which 
separates the two half-planes. 


2.5. Augmentation rules: maintaining w = 1 

It should be noted that the condition of continuity in subscripts between consecutive Ts 
(as set out in section 2.4) is not sufficient to ensure that any polymer loop described by 
a string of Ts with this unbroken property need be topologically equivalent to the basic 
loop with w = 1. This is easily seen from the string (TcTc)"', tc > 1, which is clearly 
unbroken, but is not topologically equivalent to the case w = 1. Indeed, only strings 
that are derived from each other in very specihc ways represent the same topology. To 
illustrate this, we now summarise some elementary inferences derived from the examples 
above: 


(i) A closed loop with ta = 1 in its simplest form (i.e., with the minimal number of 
piercings of the = 0 plane) is represented by the sequence 

TcTc. ( 10 ) 

(ii) A closed loop with tc > 1 in its simplest form (i.e., with the minimal number of 
piercings of the z = 0 plane) is represented by the sequence 

(T,T,)-. (11) 

(iii) The basic sequence T^T^ can be augmented (“dressed”) as in Figure [^according to 
the lengthening rule 

T, ^ TeT,T,. (12) 

(iv) A further augmentation procedure, as shown in Figure is described by the 
replacement rule 

T, ^ T,T,Te. (13) 

It is clear that compound rules arise, namely 

T, ^ T,T,T, (14) 

and 


T, ^ (Te)”T,(T,)T (15) 

We note here that ( |14D and (15) essentially encode group relations of the 
braid group B 2 — see [Appendix A) Any other alteration of the polymer strand 
through insertion / alteration of Ts that is not of type (12) or (13) (or equivalently (14) 
or (pffij)) would necessarily either break the strand (see previous section) or increase the 
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winding nnmber (see (0), This implies that the two rules ( |I4| and ( [l5| ) above capture 
all possible ways of generating loops that are topologically equivalent to the basic loop 
shown in Figure As in sections 2T and 2^, calculation of the partition function 


requires integration over various degrees of freedom. 


2.6. What sequences are valid for w = 1? 


Valid sequences generated from the simplest form TfTc {w 
must have the following properties: 


1) according to section 2.5 


(i) in order for the loop to be closed, the hrst and last index must be equal (where 
cyclic permutations of sequences are equivalent) - see equation ([^ as an example, 

(ii) for the same reason, the second index of any T in the sequence must equal the hrst 
index of the next T, 

(iii) the total number of Ts in the sequence must be even (since the basic undecorated 
closed loop has two terms, and both augmentation rules keep the total number 
even), 

(iv) the number of TcS must be even, 

(v) the number of TgS must be even, 

(vi) the string must be algorithmically reducible (this is dehned in the next section). 


2.1. Algorithmic reducibility of valid strings for w = 1 

Let us dehne the sets of generic functional units / substrings counting either even (G) 
or odd {U) sequences of T^: 

Gn = T,{Tgf^, n = 0,l,2 ,... (16) 

and 

= T,(r,)2-+\ n = 0,l,2,.... (17) 

Any string of Ts could now be rewritten as a string of Gs and Us. For the basic unit for 
w = 1 we may write TcTc = GqGq, with another example being TcTcTsTsTgTcTgTcTsTs = 
GqUiUqGi. Since the substring TgTg may be trivially inserted or removed in any 
sequence (see 0) , it is clear that 

Gn -H- Go (18) 

and 


Un ^ Uq. 


We further infer from (14) and (15) that 


A Y ^XY 


(n = 1, 2,...) V substrings X, Y 


and 


G^GoG^oGo 


(19) 

( 20 ) 


(m= 1,2,...). 


( 21 ) 
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We define a given string to be algorithmically redncible if the following procednre is 
possible: 


(i) apply (18) and (19) to simplify the string wherever possible, 


apply (20) to simplify the string wherever possible, 


(hi) now apply (21) to simplify the string wherever possible, 
(iv) repeat nntil only the fnnctional nnit GqGq remains. 


Strings that are algorithmically redncible in this manner are topologically eqnivalent to 
the basic nnit TcT^ which represents w = 1. For w > 1 the string (Go)^^ wonld remain 
in step (iv) after complete application of this procednre. 


3. Partition function 

The fnll partition fnnction for a given winding number is now given by the integrals 
over the sums of all the conhgurations that are compatible with the winding number. 


3.1. Summing over diagrams 


The rules by which moves are produced do lead to all possible configurations permissible 
as described in section]^ The corresponding sequence of Ts represents the statistical 
weight for each configuration. In order to enumerate the valid sequences correctly, 
each distinct configuration needs to occur exactly once in the partition function. 
(Alternatively one needs to be able to determine the correct multiplicity for the crossings 
in order to sum the appropriate terms in the partition function correctly.) 

For completeness we state here once more the rules from equations (14) and (15), 


(i) -)■ T^TsTs and 

(ii) Tg —)■ TcTgTc. 


The first rule adds loops of the type Tg in even multiples and the second rule is 
responsible for the addition of new terms in T^. It is simple to see that different sequences 
of applying the rules (i) and (ii) above, on different elements, can lead to configurations 
that are identical. This has obvious implications in writing expressions for the sum in 
the partition function. Here we investigate whether a scheme by which enumeration 
or an approximate enumeration are possible. (The explicit procedure can be compared 
to the configurations produced by variations of the rules and checked for repeats using 
simple algorithms in Mathematica.) 

As already explained in the previous section, the first basic consequence of the rule 
(i) above is that any even(odd) sequence of same-side crossing terms Tg can be extended 
repeatedly by a double Tg to an arbitrary degree. In this sense it is possible to use a 
compact notation for any sequence of terms in Tg and T^ by a prescription that indicated 
whether any two consecutive TcS are separated by an even or an odd number of Tg terms. 
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We utilise a notation that writes either no or one and implies the extension of the 
rule (i) summation by eventually including the factor 

1 + T,T, + T,T,T,T, + ... = (1 - T,T,)-'. (22) 

In this sense the application of rule (i) is almost trivial except when it is combined 
with rule (ii). One can hence go ahead to introduce new terms by including all the 
possibilities for odd or even expansions of Tg and then complete the series above after 
all other conhguration types have been introduced. 

It is instructive to write down a hybrid composite of rules (i) and (ii): 


(i’a) Tc —>■ TcTgTcTgTc 

(i’b) Tc TcTcTgTcTg 
(ii’) Tg ^ TcTgTc. 


We note here that the two parts rule (i’) can be interpreted in two ways: either the 
sequence TgTfTgTc is appended to the right of the original crossing or, the sequence 
TfTgTcTg is added on the left of the original Tc- We choose the the hrst of these two 
conventions since rules (i’a) and (i’b) produce equivalent conhgurations under the cyclic 
property - see [Appendix B 

Consequently, the basic winding number expression can be expanded without 
repeating conhgurations under rule (i’a) 

Tc (1 + TfTgTcTg + TfTgTcTgTcTgTcTg + ...) x ... 

= (Tc (1 - TcTgTcTg)-^ Tc (1 - TcTgTcTg)-y . (23) 


However, this clearly does not represent a sum over all possible conhgurations, since rule 
(ii’) has not been completely applied. In principle, all conhgurations should be given by 
repeated applications of the rules to all newly introduced parts of terms. The partition 
function using rule (i’a) as depicted above clearly does not repeat any conhgurations, yet 
does not produce all permissible conhgurations. We use this to calculate an approximate 
partition function 


y(w) 

■^appxl 


Tc 


2w 


(1 - T.T.) [1 - t.t, (1 - r,r,)-‘ T.T, (1 - T,T.)-'] 


(24) 


Integration over relevant degrees of freedom of this expression is implied. Since the 
weight of each conhguration in eq. (24) is the same as in the complete sum for the 
partition function the complete partition function for winding number w given by 
is related to the approximation as follows, 

42x1 < Z'”’- (28) 


(We note that careful implementation of rule (ii’) on a subset of Tg terms above will 
lead to an ever better lower bound than .^appxiO 
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Another interpretation of iterative application of rules 
definition of two coupled effective terms 




-1 


“’eff 


t,t; 


effyefF 


'^611 _ n~^i 

s — s ~r -L c 


efFyefFyeff 


and (ii’) is given by the 

(26) 

(27) 


Here explicit evaluation shows that the scheme eventually does lead to repetition of 


some conhgurations, but all conhgurations are produced when combined with (22) at 


the last step. The partition function calculated using the recipe in (26-27) is dehned by 

_ j’^efFrynefFl 


7M 

'^appx2 


ny]' 


(28) 


Now since eq. (23) leads to a partition function .^appxi with correctly weighted, yet 


fewer conhgurations, and equations (|26|-|27|) yield a partition function with all 


yet some multiply occurring conhgurations we know how the true partition function is 
bounded 


yM ^ yM 

^appxl — ^ — ^appx2- 


(29) 


In principle these two approximations are calculable in the scenario of a polymer loop 
winding around certain obstacles in the plane, as described in sections and and 
can be used to understand upper and lower bounds for free energy associated with 
a particular winding number. We calculate only the lower bound .^appxi here, as the 
nonlinear coupled equations (26) and (27) pose formidable challenges. As stated, the 


approximations still need to be integrated over the relevant degrees of freedom, as will 
be set out in section 13.41 


3.2. Counting the number of crossing or same-side terms 


In either of the suggested approximations for the partition function (see (24) and (28)) 
it is possible to include generating terms that may be used to calculate the number of 
Tg or Tc terms. If in these summations we simply replace Tg —)■ e^'^Tg and Tc —>■ e^^^Tc, 
then we may calculate the average number of crossing terms as 




Aiogz(”i(j„j,) 

dgc 


-I 9c=9s=0 

and the average number of same-side terms as 


{Ns) = 


■^logZ^^\gc,gg) 
[9 9 s 


(30) 


(31) 


9c=gs=0 


3.3. Probability distribution of a flexible polymer in half-space 

We still have to assign a statistical weight to each string in the summation over all 
diagrams. We proceed to do this for a flexible polymer. In principle other polymer 
variants could be described by the formalism up to this point, but the form of the 
probability distribution would be different. 




















Conservation of polymer winding states: a combinatoric approach 


14 


For a flexible polymer we treat the sub-arcs (labelled by the various Ts) as 
random walks confined to half-spaces. In polymeric systems with suitable solubility 
and flexibility conditions this is, of course, a reasonable assumption (2^. The notion of 


confined random walks is certainly not a new one. In a 1943 review Chandrasekhar 26 


pointed out how one-dimensional random walks with reflecting and absorbing boundary 
conditions may be treated. Naturally boundaries change the probability distribution 
for random walks. A reflection off such a boundary implies that a walker must 
necessarily retrace its last step. An absorbing boundary, in contrast, would prevent any 
further displacements. Consequently the probability distribution of a walker confined 
by a reflecting boundary is obtained by adding an “image distribution” to that of an 
unconhned walker. This accounts for additional possible paths to a given end point, 
stemming from the reflecting boundary (these paths may be viewed as mirrored paths 
in the excluded region). On the other hand, the distribution of a walker confined by an 
absorbing boundary is obtained by subtracting a similar mirror distribution. This, in 
turn, accounts for the exclusion of trajectories that terminate on the absorbing boundary. 

This discussion may be extended to a random walker restricted to a half-space in 
three dimensions (see, for instance, the article of Slutsky 27 where a similar “method 
of images” is used). We shall draw on these notions in order to assign appropriate 
statistical weights to the sub-arcs mentioned in the previous section. To this end we 
make the following assumptions (as illustrated in Figure]^: 


(i) there exists a finite minimal length scale (such as a bond length or Kuhn length) 
in this polymer system, 

(ii) for every (sub)sequence T^yTy^ (for any x,y, z & —}), there exists a trans-plane 

polymer segment of length 2e that is normal to the plane and connects the two 
sub-arcs between the > 0 or the z < 0 half-spaces, 

(hi) thus any given sub-arc begins and ends at a distance e from the xy-plane (see 

(iv) each polymer sub-arc is modelled as a random walk constrained by an absorbing 
boundary plane to either the z > 0 or the z < 0 half-spaces, 

(v) such a random walk begins and ends at rj) = (xo,2 /o,?7 e) and r = {x,y,rie) 
respectively (here p = +1 for the z > 0 half-space or p = —1 for z < 0), 

(vi) the random walks are fully flexible, and each has a variable arc-length Sj. 



The assumption of an absorbing boundary is based on the fact that we are interested 
in the two piercings that a sub-arc makes with the plane, since it is there that one 
particular sub-arc ends and another one begins. In reality, the trans-plane connecting 
segment should be free to take on any orientation. Our approximation that it is normal 
to the plane should be a small correction for a sufficiently long arc segment. It is clear 
from Figure that, barring constraints due to the total length of the polymer, the y 
components of the beginning and end vectors of a random walk could take any value. 
The X components, however, are restricted either to the interval (0, oo) or to the interval 
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Figure 6: Sub-arcs as random walks that begin and end at a distance e from the plane. 


(0, —C)o), depending on whether they fall on the -|- or — side of the rod. Naturally this 
will constrain the integration bounds for x in the partition function accordingly. 

Therefore, for a Gaussian chain with all lengths expressed in terms of the Kuhn 
length, we may now assign the corresponding probability distribution to a typical sub¬ 
arc as set out in (2^, 


r(f,ro,s,e) = (27rs) 
= (27rs 


^-(r-ro)^/2s _ -(r+n))2/2s 


-3/2g-[(x-a:o)2+(j/-yo) ]/2s 


-eCs 


(32) 


where s is the arc length of the polymer sub-arc between vq and r. The 2 ; dependence 
vanishes due to assumptions and ( pfi] ) above. Here we have interchanged the discrete 
number of steps N for the arc-length variable s through appropriate re-scaling. Finally, 
given a polymeric system with a sufficiently small Kuhn length, we may assume e <C 1 


and Taylor expand the exponential in Equation (32) to obtain 

pi 


T 

-L n 


r,ro,s,ej = 


a/ (27r)^s^ 


^-[{x-xoY+(y-yor]/ls 


(33) 


The subscript p here simply refers 


consideration, as described in section 2.4 


to the “parity” of the particular T under 
This merely indicates whether the x co¬ 


ordinates are on the same side of the rod or not, which will determine the integration 
domains for the x co-ordinates in the partition function. 


3 . 4 . Partition function for w = 1 


Using section 2.5 we constructed a symbolic summation of all possible configurations 


for tn = 1 in two possible approximations, (24) and (28). We now proceed to use this 


summation in order to write the partition function for this system. Let us denote the 
set of all valid conhgurations as A. Supposing that the total length of the polymer loop 
is L, we note that 




X6A 


( 34 ) 
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i.e., the total partition function is simply the sum of the partition functions for all valid 
conhgurations for w = 1. A typical valid conhguration y which has N piercings of 
the a;|/-plane is represented by some sequence (of length N) of Ts that adheres to the 
conditions in Section |2l6 , and simply corresponds to one of the terms in the summation. 
Let us make the notation in equations ([^ and (|^ more concrete: the partition function 
corresponding such a generic sequence would be 


= / dX dY dS 6ixo - XN)S{yo - 2/iv) 
J 0 J—oo J e 

N N 

'^Pj (^i—1) Vj—^ Vji Lf 

j=l k=l 

The condensed notation implies 

dX = dxo dxi... dxN, 
dY = dyo dyi... dy^, 
dS = dsi ds2 .. .dsN, 


(35) 


(36) 


and the functions Tp in the integrand are each of the form (33). Naturally there is one 


less s integral than for x or y, since one arc-segment connects two planar points. The 
first two delta functions ensure that the strand is closed: the first and last x and y 
co-ordinates must be equal. The length Si of each of the N sub-arcs is bounded from 
below by the minimal length e - hence the integration bounds on the s integrals. It is, 
however, necessary that these lengths add up to the total length L of the entire polymer. 


This constraint is enforced by the third delta function. As set out in Section 3.3, the y 
co-ordinates of each strand are not constrained. For this reason they are integrated over 
the whole axis. From ( [3^ it is clear that the Ts are symmetric under the exchanges 
(xj_i — Xi) —)■ —(xj_i — Xi) and {yi-i — yi) —)■ —(l/i-i — yt)- Since the x co-ordinates 


are constrained to one half of their axis (as set out after Equation (33)), we need to 


distinguish between the terms that cross over the rod iTf) and those that remain on 
the same side of the rod (T^). For T^, the two x argument have the same sign and the 
function depends on ±(xj_i — Xj). For Tc, the two x arguments have opposite signs and 
the function depends on ±(xj_i -|- x,). We may thus change all x integrals to run over 
(0, cxd) under the condition that 

^2 


Ts{xi-i, Xi) = T{xi-i - Xi) = 


Tc(xi_i, Xi) = T(xi_i -F Xi) = 


V(27r)3sf 




^-[{xi-i+XiY + {yi-yi-lY]C^i 


^/W) 


3 


(37) 


i.e., with these integration bounds the same-side contributions depend on the difference 
between their x co-ordinates, whereas the crossing contributions depend on the sum 


(compare to (33)). For the reasons set out above, it is clear that translational invariance 
holds for the y co-ordinates but not for the x co-ordinates. In [Appendix C| we outline 


how some integrals in (35) may be diagonalised using a Laplace transformation in the 
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length of the polymer and Fourier transformations of the y co-ordinates. The result is 
the Laplace transformed partition function for some configuration y, 

/»oo /•no ^ 


zy '*(«)= / dXg(xa-XK) / dk 


(38) 


i=l 


where the superscript implies that the Ts from (37) have been Laplace 

transformed and Fourier transformed in y. We have thus obtained a diagonalisation 
for the s and y co-ordinates. What remains are the integrals over the positive (real) 
X axis, and the integral over the Fourier variable k. We shall deal with these integrals 
for two cases. First we consider restricting the x integrals to a narrow slit, thereby 
effectively constraining the polymer to be wound through two slits in the plane. In this 
case no x integration is necessary, and only the k integral remains. Secondly we shall 
outline possible approximation schemes to deal with the general case. 

As stated, the complete partition function is found by the summation over 

various diagrams. This summation is approximated by (24) or (28). The integrand 
in this partition function may then be repeated w times to obtain the approximated 
partition function for higher winding numbers w > 1. 


4. Specific case: polymer wound through two slits 


Let us consider a polymer looping around two slits, the inner edges of the slits separated 
by the distance d and each slit with a width A, as shown in Figure [T} W e note that 
similar scenarios of confined wound polymers have been considered in 


19 



Figure 7: Constraining the polymer to two narrow slits in the plane. 


We can compute the average number of the types of arc elements in this scenario 
using the formalism developed in section 3A By winding the polymer around a double 
slit geometry, the combinatorics clearly remain unaltered. The partition function for the 
chain now has integration restricted over the domain Xi G [d/2, d/2 + A]. This means 
that equation (38) takes the integral form 

l+A 


N 


^\t,A,d)= / dX 6 {xo - xn) / dk Y\_T^’^{xi_i,Xi-,k-,t).{39) 


2=1 


We begin by considering the case of zero slit-width. 
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In the limiting case where A —)■ 0 such that each T has exactly the same x-value 
Xi = (i/2, Vi G {0 ,...,A}, the partition function becomes especially simple as the 


operator sums of eq. (24) or (28) now become simple algebraic sums. (This is the scenario 
where no integration is necessary, and conformations are simply summed because the 
problem is “diagonal” already in the simplest terms. The next section will deal with the 
extension of this to narrow and easily integrable slits.) Choosing the hrst approximation 


(24) we obtain 


}{w=l) 

■^appxl 


(f, d) = 


dk 




[2 + 


^ 2 ’ 


(40) 


since the operations of summing over various configurations and integration over k may 
be exchanged. 


Turning to equation (37) we see that Fourier transformation in y leads us to the 
following two cases, 

t.2 „ 

rpF ( 1 1 \ 

T, [d; k; Si) = 


2s 


(41) 


and 


Tc{d;k;si) = 






2 s 


(42) 


Since e is non-zero and finite, these expressions are well-defined. We now need to perform 
the Laplace transforms of each of these. 


TY^{d]k]t)= / dsi e 


-.n \/3e2e' 


3^2 


2 s? 




-|e(fc 2 + 6 t) 1 / 1 

-/-g + Qt)r U,-{k^ + 6t)e 


and 


TY^id;k;t)^ / dsi e 

JdP 


f— r, 


2s 


(43) 


(44) 


In (43) the answer contains an incomplete gamma function. The integral in (44) has an 


approximated lower bound of d? (recall that s* and d are dimensionless). In principle 
this bound should be e (i.e., of the order of the Kuhn length). If the polymer were 
inextensible, the minimum arc-length for should be d. Although we deal with a 
Gaussian chain here, this approximation is reasonable since the integrand is dominated 
by Si> d^. We approximate 


{d]k]t)’Fs-e / dsi e 

2 Jd? 


-.4n\/3e‘e 


2 „-- 


2 s 


( 45 ) 
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where s* = ds se 6‘')/(J^ds e s'*) = 


-V'l — 6+d^k^ 


is a constant valne that captures 


'd^ / ]^2 

some of the k scaling of the answer. Numerical verification shows this approximation 


to perform very well for various ranges of t and k. The integral in (45) may now be 
evaluated explicitly, 


T^^’^{d-k-t) = 


- 


2 ( 6 +d^fc^) 


X 




- + 6t) r 1 ^ 0 , {k‘^ + 6t) 


.(46) 


We may now insert the answers (43) and (46) into (40) to obtain two approximations 
for the partition function. Naturally the integrand above is some very complicated 
function of k. As is verifiable numerically, however, a saddle point approximation is 
reliable for various ranges of t and d. We omit the cumbersome explicit form of the 
result. 


The average length of the loop for w = 1 may be calculated from (40), 


(r)W = -|iog 


}r{w = l) 
■^appxl 


{t,d) 


(47) 


In Figurej^we see that {L)(t) is a concave function for various values of the slit separation 
d. Here we have set e = 1 for convenience; this convention is used henceforth. It is also 
clear that small Laplace parameters and longer length scales correspond, particularly 
for small d. As we increase the size of d, we note that {L) seems to strive asymptotically 
to increasingly large values. This is simply a manifestation of the fact that the non-zero 
slit-separation implies a minimal length scale for the polymer loop. 

We thus have a “lookup table” that allows us to associate an average length of the 
polymer to a particular Laplace parameter (of course this is not the inverse Laplace 
transform, as would ideally be the case). 

im) 

80 h 


60 


40 


20 


Figure 8: Average length of the loop as function of the Laplace parameter t, calculated 
according to equation (47). Parameters: d = 1 (solid), d = 3 (dashed), d = 5 
(dashdotted). 
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It is further interesting to ask what is the relative weight of the undressed term / 
basic loop TfTc (see Figure in the summation over all valid diagrams for ta = 1. To 
this end we simply look at the probability for this configuration, 

iz.dk 

zltCit.d) 

We see in Figure that, for various values of d, this probability is unity for snfficiently 
large Laplace parameters. Throngh Figure!^ we may thus identify length-scales (for 
various slit separations) at which the basic (nndressed) loop provides the dominant 
contribntion to the partition fnnction. This makes sense physically, since large Laplace 
parameters correspond to short length-scales. Naturally the afore-mentioned minimal 
length scale set by the slit separation implies that as soon as the length of the loop 
becomes small enongh, the basic configuration consisting of two crossing terms will be 
the dominant configuration in the partition fnnction. 



P{T^ 


P(Tc^) 



Figure 9: Probability for the configuration as a fnnction of the Laplace parameter 
t (calculated according to eqnation (48)). Parameters: d = 1 (solid), d = 3 (dashed), 
d = 10 (dashdotted). 


The average number of crossing terms and same-side terms may also be calculated 
as a function of the Laplace parameter according to eqnations (30) and (31). From 
Figure 10 we note that as soon as the Laplace parameter is snfficiently large (i.e., the 
polymer is typically short), there are exactly two Tc terms and zero Tg terms present. 
This agrees with the previons conclusion: at short polymer lengths, the nndressed basic 
term TcTc dominates the snmmation over valid diagrams. We note that as the slit 
separation d is increased, this nndressed configuration dominates at decreasing t, i.e., at 
longer polymer lengths. This can be related to Figure where we observe the minimal 
length of the polymer increasing with increases in slit separation. 

As is to be expected for snfficiently small polymer length-scales, the results are not 
particularly sensitive to which approximation is used to approximate the summation. 
Indeed, for short or stiff polymers it shonld be snfficient to generate the first few terms 
(valid seqnences) explicitly as an approximation to the complete partition fnnction. 
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m(t) 



(b) Average number of same-side terms 


t 


Figure 10: The average number of crossing and same-side terms as functions of Laplace 
parameter, calculated numerically from (30) and (31). Parameters: d = 1 (solid), d = 2 
(dashed), d = 3 (dashdotted). 


Lastly we consider the free energy of the system for various values of slit separation. 
This quantity is simply the negative logarithm of the partition function, and is plotted 
parametrically as a function of the average polymer length in Figure It is clear 


that as the polymer is made shorter and approaches the minimal length-scale set by 
the slit separation, the free energy increases sharply which is compatible with the sharp 
decrease of entropy experienced by the Gaussian chain. 


F 



Figure 11: Parametric plot of free energy dependence on average polymer length. 
Parameters: d = 1 (solid), d = 2 (dashed), d = 3 (dashdotted). 


4-2. Finite slit width: A 7 ^ 0 

We shall briefly outline a possible approach to solving the case for non-zero slit width 
in (39). As a hrst order approach it would be sensible to decouple consecutive Ts 

— +A 

completely and to replace the various x integrals (ff dxi) with a localisation 
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approximation, 


and 


dx 


dx 


' —oo 


A 


dx' ^ 
A 


exp 


exp 



for Tg and Tc respectively. The analogues of equations (|41|) and (|42|) now become 

Tf{d;k;si) = 


y/Se^e 


2 „-- 


2 si 


Si + 2A 


and 


T^{d-k-sf) = 


^/3r2 


6 e 


■ 2s,- 


2 s 


+ 2A 


Zd^ 

'2(2A^+Si) 


(49) 


(60) 


(51) 


(52) 


Naturally these expressions reduce to the case of zero slit-width if A —)■ 0. Again we 
require Laplace transformations in order to use (24), and we must do the k integral as for 
(40). This may be approximated in various ways. A Taylor expansion of (51) and (52) 
to 0{A), for instance, allows us to repeat the analysis from section 4.1 without many 


modihcations, yielding an approxima 


(We may recover (1^ through limA^o d, A) = Zf^pAi A, d).) 


ion for Laplace transformed partition function. 


{w=l). 




One may now calculate, for instance, the force exerted by the slit as the derivative of 
the free energy. 


d 


f{t,d,A) = (^-log ^zipJ^\t,d,A) 

With the aid of Mathematica we may “invert” the Laplace transformation numerically 
to obtain the Laplace parameter as a function of the average length, t = t{{L),d,A). 
This allows us to plot, for instance, the force as a function of slit separation for a hxed 


(53) 


(L), see Figure 12 





Figure 12: Force exerted by the slit as a function of slit separation d, for A = 0 and 
a hxed (L) = 20. The ratio exhibits a peak corresponding to the minimum of the 
force. Compare to Figure |8| to see why d > 3.1 is excluded. 
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Despite the somewhat crude simplifying assumptions made (using independent 
localisation of each arc), expected physical aspects are well captured. For a fixed average 
polymer length, there is a sign change in the force at some slit separation. This makes 
sense: for sufficiently small d the slit has a “compressing” effect on the polymer, and for 
sufficiently large d the polymer is “stre tched ”. Again the competition of the two length 
scales d and (L) is manifest. In Figure 


12 b 


we plot the ratio as defined in (30) and 
(31) for a fixed (L). Consider the trend in this plot as we decrease the slit-separation d. 
For large d all of the polymer is in the crossing terms and (Ng) ~ 0. As we decrease d the 
fraction of same-side terms increases: more of the polymer’s length is free to occupy the 
slits. As d is decreased even further, the ratio begins decreasing again: at some stage 
sufficiently much polymer length is free so that additional crossing terms may arise, 
thereby decreasing (Ng) and increasing (Nf). Corresponding behaviour of the force is 
evident in Figure 12a (also compare to Figure [l0|. 


5. General case: outline of solution strategy 


Suppose we return to the general form of the partition function in equation (38). For 


the general case of an infinitesimally thin rod in the plane, the remaining integrals over 
X are not as easy to diagonalize. Instead of such a rod we shall consider a flat slab of 
width d lying along the y axis in the xy plane; see Figure This amounts to modifying 
the two-slit scenario of previous section by removing the “outer barrier” of each slit, 
thereby changing the x integration domains to (|, cxd). The limit d —)■ 0 represents the 
original scenario of Edwards’ rod in the plane. 

Recall that the terms in the partition function consist of multiples of Tc with even 
or odd geometric series of Tg, i.e., Tc(l — Tf)~^ or TcTs{l — We approximate 

these contributions to the partition function in two steps: 

• address the sub-sequences of single-side contributions that originate from 


augmentations of the type (12), i.e., 

Tc{x + Xo) Tg{xo - Xi) ... Tg{xm-i - Xm) Tc{xm + x'), and then 

'-V-' 

m single—side terms 

approximate the crossing terms T^. 


5.1. Approximation of Tg sequences 

Let us begin by considering a sub-sequence of m single-side contributions in the integral 
(38) between xq and x^, say. For notational convenience we omit k and t dependence 


(which is the same in all terms), but recall that the k integral still remains. Let us 
consider the case where we have not yet performed the Laplace transformation w.r.t. 
L, and define 

poo poo 

Tfg{xo,X^]t,{Si}) = / dxi... dxm-i 


'd/2 


ld/2 


Tf (xo - Xi, k] Si)... T[{x^_i - X^] k] s^). (54) 
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Here each has been Fourier transformed in its j/s but not yet Laplace transformed, 
i.e., 


T^(t- T- 1-h- s i — 

s J — 2 ^ ’ 


( 66 ) 


as is easily seen from (33). This implies that we may write 

m—1 


^ m—1 


Y[Tf{xi,Xi+i]k]Si) = JJ 


■hf _l.‘2 


.3/2 


—k^Si/d 


m 1 ^ 3(xj-x^^lP 


2 = 1 


. 2=1 


n, 

.i=i */ 


1/2 


e 


(56) 


We note that second product is simply one of Green functions for a one-dimensional 
polymer chain, 


C (r- T- 1 - s i — -1,(xi-xi+if/2s, 

^'^i+li ^i) 1 / 2 ^ 


(57) 


The restriction imposed by the integration bound on each intermediate Xi (where the 
polymer pierces the plane) implies that none of the piercings may enter the excluded 
region of the bar in the region x G (—|, |)- This however does not preclude any other 
part of the polymer arc to cross over this region, as illustrated in Figure [T^ 


d d 



Figure 13: Piercings are excluded from a slab region on the y axis, but the remainder 
of a polymer arc could still cross over this region. 


We proceed with our approximation by assuming that indeed none of the segments 
in an arc may cross the excluded region, i.e., all x co-ordinates of an arc (and not just 
those of its piercings) must lie in the region (|, cx)). Naturally this approximation would 
exclude several possible conhgurations, and we are signihcantly under-estimating the 
true entropy of the polymer arc. As such this would provide a sensible upper bound 
for the free energy of the system. This case should be adaptable to the scenario of a 
polymer in a cavity from 19 through a similar argument in a cylindrically symmetric 
setting. 

Under this assumption the product G^ixi, Xj+i; Sj) reduces to a single chain of 

length that nowhere crosses into the forbidden x domain, i.e., the polymer arc 

is effectively restricted to a quadrant around the forbidden region. For such a restricted 
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random walk we may once more use Chandrasekhar’s argument of mirror images to 
construct an effective Green function, 

m—1 m—1 


2 = 1 


2 = 1 


d.. 






d 


m— 1 


2=1 


- -) + (i., 


d 


m—1 


( 68 ) 


2=1 


Here the superscripts on the Green functions refer to the respective x integration 
domains. We have thus replaced the sequence of G^s with an effective Green function 
by making use of the Markov property of random walks after having extended the x 
integration to all space by making use of Chandrasekhar’s argument. This is manifest 
in the last line above of equation (58) where the Green function with the reflected 
end co-ordinate has been subtracted. This means that all intermediate x integrals 
have vanished, and only the integrals over the beginning and end x co-ordinates of 
G'^x~d^‘^\xQ^Xm'iYlT=i^'i) remain. In this approximation the dependence on m has 
completely disappeared out of one part of the expression in equation (56) yielding 

m—1 


' m—1 


M 


m—1 




2 = 1 


,*=i s 


2 = 1 


What remains is to perform the integrals over the s co-ordinates, with the relevant 
cut-off on the lower integration bound, 

/•OO /•OO 


Tf^ {xQ,XmG) = / dxo dXr, 


n 




m—1 


Af 


m—1 


n x„; Y1 »<)■ (60) 

i=l i=l 

This effective quantity depends only on the first and last co-ordinates, x and x', and on 
the number m of same-side steps in the sub-sequence. A summation over m will lead 
to a geometric series x') that may now be included between Tc terms in a 

chosen approximation scheme for the summation over valid sequences. 


5.2. Approximating the Tc terms 

The second step is to deal with rod-crossing terms of the form Tc{x,x'). These terms 
are responsible for the localisation of the polymer around the rod at these points. We 
propose an approximation that decouples x and x', 

Tc{x,x) ^ f{x)f{x). (61) 

This could be done in several ways, and would lead to a complete diagnoalisation of 
the integration. For a short polymer one may assume that the crossing terms will be 
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localised close to the boundary of the slab, i.e., {x) ~ For other length scale regimes 
this average localisation could also be guessed as a function of the Laplace parameter t 
or solved in some self-consistent manner. 

A slightly more general approach would be to assume that the piercings of the 
crossing terms are localised by Gaussians around x,x' = ^. Naturally one would expect 
a Tc term to show a corresponding decline if a; G (|,cxd) and x' G (f,cxD) are moved 
further from the origin, since it is a Gaussian in x + x'. The Laplace parameter t is 
related to the inverse of the arc-length of the crossing segment, and should be indicative 
of this localization. A sensible approximation for the x dependent part would thus be 
to say that Tc localizes x and x' independently (i.e., we decouple the function into two 
independent Gaussians), 




[x 


x') 


UqC 


-fG-fd-fG'-f) 


(62) 


where the localization parameter q is treated as a guess to the localisation length scale. 
Here Mq is the normalisation such that 


dx 


dx' Ar„e-f = 1 


(63) 


One could also determine the strength of the localisation using, for instance, a 
variational calculation. The result could now be combined with those of section IaD and 
section 3.1 in some approximation for the partition function. 


6. Conclusions and closing remarks 

The problem of winding a polymer around an inhnitely long obstacle was addressed 
by labelling conhgurations according to sequences of sub-arcs constrained to a half¬ 
space. By considering arcs that cross the obstacle or remain on the same side thereof, 
combinatoric rules were derived that allow winding number conserving augmentations 
of these sequences. These augmentations rely on topology conservation through type 2 
Reidemeister moves. Properties of valid sequences were identihed and an algorithm was 
presented for discriminating whether a given sequence is valid. 

Two possible approximations for the partition function were found by considering 
summations over valid conhgurations, one bounding the true partition function from 
above and the other from below. Given a particular choice of polymer variant (Gaussian, 
semihexible etc.) a statistical weight may be attached to each sub-arc. Through a 
series of diagonalising transformations the lower bound approximation for the partition 
function was written as an integral of these statistical weights. 

For the case of a Gaussian chain a specihc statistical weight was assigned, and the 
partition function for w = 1 was approximated for the case of windings through two slits 
in a plane. For zero slit-width it was found that the basic undressed loop dominates the 
partition function when the polymer is short, as expected. Various expectation values 
were calculated and the free energy was plotted for different slit separations. The results 
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make good physical sense for various average polymer length scales, lending credibility 
to our approximations. 

The case of non-zero slit width was treated in a localising approximation. The force 
of the slits on the polymer shows an expected sign change as slit separation is varied 
for a hxed average polymer length. Correspondingly different length scale regimes arise 
that determine what fraction of the polymer is in crossing terms or in same-side terms. 

Only the case w = 1 was treated in detail. Since the partition function arises from 
summing integrals of products of functions, the partition function for higher winding 
numbers may in principle be constructed from higher order basic loops by our rules. 
This implies that some statistical quantities calculated from logarithms of -^appxi (e.g., 
{L){t), (Ng), (Nc)...) scale with the winding number of the system. In our lower 
bound approximation scheme, at least, higher windings entail repeats of the integrand 
a corresponding number of times. The upper bound approximation scheme could also 
be useful in some regimes since the hrst few terms in the summation over conhgurations 
dominate the partition function for certain length scales. 

Lastly, suggested approximations were sketched for a slab-like obstacle in the plane. 
It was observed that the latter case should be relatable to other conhned geometries 
studied in existing literature. Our approach allows different and intuitive calculation 
strategies for partition functions of wound polymers. 
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Appendix A. Comments on enumeration and braid groups 


We explain here briefly the relation of our enumeration procedure of section to some 
properties of braid groups. Since we wish to calculate partition functions and other 
averages for polymers of a given winding number we illustrate why braid group properties 
alone do not suffice to couple polymer degrees of freedom to a particular conhguration. 

A braid group Bn describes braids of n strings, and has n — 1 group generators 
{cTi, (72,..., Cn-i}- The generators obey the relations 


^ 1, 2, . . . , 77 . 1, 

^j 1 ^ j\ ^ 2 , 

CTiCTj”^ = cr“Vi = e. (A.l) 


Here e is the identity element. One may construct words from this set of generators 
and their inverses; these words correspond to braids of the n strands. (See, for 
instance, 20-23 for further details.) For any given word the application of the braid 
group relations (A.l) may be used to obtain a minimal (irreducible) form of word. 
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The scenario of winding a strand aronnd a rod could, of course, be viewed in terms 
of the braid group for two strands, B 2 . This group is “trivial” in that it only has one 
generator a, rendering the hrst two properties of (A.l) irrrelevant. In Figure Al we 



Figure Al: Braids of two strands in analogy to winding scenarios from section 


illustrate the B 2 words analogous to some sequences of Ts from section The braids 
shown there are projections of the winding scenarios along the plane in which the rod 
lies. Clearly the group properties thus encode similar conditions on reducibility as set 
out in section Since we consider words of even length, the winding number of a 
particular B 2 braid may be obtained by halving the number of a generators remaining 
in the corresponding word that has been fully reduced. For any word of even length, 
the winding number is clearly related to the difference in number of as and a“^s. 


w = 


#(^) - 


(A.2) 


(Compare to the power of the Alexander invariant in equation (2.58) of |^.) 

However, our aim is to couple polymer degrees of freedom (statistical weights) to 
these words. These weights are simply the probability distributions of connected arcs 
constrained to half-space, viz. the Ts. It is critical to identify the correct integration 
bounds on the beginning and end co-ordinates in these weights (see sections 3.3 and 


3.4) when stringing together these Ts for a given braid / winding scenario. Indeed, the 
group relation is problematic in that it alone does not, for instance, allow 

us distinguish between the sequences TgTgTg and TflsTc. Clearly, through (14) and (15) 
both of these sequences reduce to Tg. In our partition function, however, where the 
TcS and T^s do not have the same integration bounds, these two conhgurations would 
represent different statistical contributions. 

Naturally it is possible to enumerate all unique (non-primitive) words in the “braids 
only” scenario, where the winding number constraint w = 1 may analogously be viewed 
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as enumerating all words in B 2 that reduce to aa. One could do this by enumerating 
all possible words of even length, and then enforcing w = 1 through (A.2) and an 
appropriate Kronecker delta. However, the translation to arcs with position dependence 
must be made explicit to avoid the above-mentioned problem. The group relations alone 
(without position dependence) thus do not allow coupling to the polymer degrees of 
freedom. Our scheme has thus been developed to avoid the necessity of including the 
ff{(T) — constraint in terms of an additional integral for the Kronecker delta by 

enumerating the conhgurations (with position dependence) explicitly. The sequence of 
Ts does matter since it indicates the order of coupling of positional degrees of freedom. 


Appendix B. Redundancy of rule (i’b) 

To see why rule (i’b), for example, is not necessary when using rule (i’a), we consider 
how an initial knot conhguration, given by a pure winding number such as 

= (wr (B.i) 

is modihed by the application of these two moves. Inserting rule (i’) after the second 
Tc is equivalent to inserting T^T^TcT^ after the hrst T^. Further equivalence can also be 
avoided by noting that the composite rules produce products of terms of the form 

T, ^ {TfTsTcTs + T,T,T,T,)™ . 

This can simplihed by noting that the following cross product terms arise, but can be 
derived from from (ii’) with the odd/even Tg convention 

TcTgTcTg X TgTcTgTc ■<— TfTg 

jii') ^ 

then even Ts after 1st Tc 


Appendix C. Various transformations in section 3.4 


The general form of (|35|) may be simplihed by performing a Laplace transformation. 

dsi... I dsN d{xo-XN)d{yo-yN) 


Suppressing all x and y dependences, we condense (35) to 

/*oo /*oo /*oo /*OC 


dX 


dY 


N 


Tpj(si)...Tp^(sjv) 6{J2sk-L). 


(C.l) 


k=l 


Performing the Laplace transformation (indicated by the tilde below) yields 


-Lt 


Z^^=^\t) = dL e 


= dX5{xo-XN) dY S{yo-yN)T^,{t)...T^^{t),{C.2) 
Jo J —00 

i.e., by performing all integrals over s we obtain the product of the Laplace 
transformations of the N individual Ts. The superscript Ts indicate the Laplace 
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transform of (33). It is, however, understood that all s integrals are performed from the 


non-zero lower bound e. The object (C.2) may further be simplihed in terms of the y 


integrals. To this end we perform the linear and invertible co-ordinate transformation 
A?/j = yi — yi+i, i = 0,1,... iV — 1 and R = J2iLo Uk- The Jacobian is \J\ = iV -|- 1. Now 


suppressing the x and t dependences we may write (35) as 

poo poo 

/ dX6{xo-XN) dyo... 


dyN 


Tp,{yo-yi) 


J —oo J —oo 

■ Tp^iVN-i - Vn) d{yo - yN) 


dR 


= \J\ j dX6{xo-XN) dAyi... dAy^ 

Jo J —oo J —oo J —oo 

N 

k=l 

poo poo 

= Af / dX6ixo-XN) / dk T^f{k)...Tif{k). (C.3) 

Jo J —oo 

where Af = \J\ dR is a pre-factor related to the length of the rod. This pre-factor 
is, in principle, divergent, but its contribution to the free energy is additive and thus 
irrelevant. The superscripts in the hnal line above indicate Fourier transformations in y 


of the Laplace transformations of the Ts from (33). We used the Fourier representation 
of the delta function, 

^ poo 

6{J2^yk)= dfc (C.4) 


k=l 


to do all the Ay integrals. Recalling that in (C.3) we suppressed dependence on 


the Laplace parameter t and on the various x co-ordinates, we write the Laplace 
transformation of explicitly as 

poo poo ^ 

Zi^=A^t)= dX6{xo-XM) dfc (C.5) 

Jo J—oo •_! 


as in (38). 
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